o 
o 



Oh! 



X 



hep-th/0409162, CALT-68-2520, NORDITA-2004-77 

Interaction Rates in String Gas Cosmology 

Rebecca Dano^ 
Department of Physics, 
University of California, Los Angeles 
Los Angeles, CA 90095 



Andrew R. Fre\ 



O 

^ I California Institute of Technology, 452-48 

00 



Pasadena, CA 91125, USA 

' Anupam MazumdaJil 

NORDITA, Blegdamsvej -17 
^ . DK-2100, Copenhagen, Denmark 



Abstract 

We study string interaction rates in the Brandenberger-Vafa scenario, the very early universe 
cosmology of a gas of strings. This cosmology starts with the assumption that all spatial dimensions 
are compact and initially have string scale radii; some dimensions grow due to some thermal or 



^ ■ quantum fluctuation which acts as an initial expansion velocity. Based on simple arguments from 

the low energy equations of motion and string thermodynamics, we demonstrate that the interaction 
rates of strings are negligible, so the common assumption of thermal equilibrium cannot apply. We 
also present a new analysis of the cosmological evolution of strings on compact manifolds of large 
radius. Then we discuss modifications that should be considered to the usual Brandenberger-Vafa 
scenario. To confirm our simple arguments, we give a numerical calculation of the annihilation rate 
of winding strings. In calculating the rate, we also show that the quantum mechanics of strings in 
small spaces is important. 

PACS numbers: 11.25.Wx,98.80.Cq 



* Electronic address: |rdanos@physics.ucla.edu| 
^Electronic address: frey@t Tieory.caltech.edu| 
■I- Electronic address: anupa nim@nordita.dk| 



1 



I. INTRODUCTION 

In Brandenberger and Vafa (BV) proposed a seemingly very natural initial condition 
for cosmology in string theory. In BV cosmology, all nine spatial dimensions are compact 
(and toroidal in the simplest case) and initially at the string radius. The matter content of 
the universe is provided by a Hagedorn temperature gas of strings. In addition to proposing 
a very interesting initial condition and analyzing the thermodynamics of string at that point, 
however, BV argued that string theory in such a background provides a natural mechanism 
for decompactifying up to three spatial dimensions (that is, allowing three spatial dimensions 
to become macroscopic). The BV mechanism works because winding strings provide a 
negative pressure, which causes contraction of the scale factor, as was shown explicityly in 







BV then gives a classical argument that long winding strings can only cross each 
other in three or fewer large spatial dimensions. Therefore, since winding strings freeze out 
quickly in four or more large spatial dimensions, the winding strings would cause recoUapse 
of those large dimensions. 

Ill has inspired a broad literature. One important generalization has been including 

nnriQ 

branes in the early universe gas of strings ULISLISI) IH) other spacetime topologies have 




also been considered P, 1^. Starting with a number of authors have examined the 

cosmological equations of motion app ropriate to string g well as brane gases 

fl Q H Q H 0, Q H Q 21, 22, in particular, 0, Q 



showed that 



interesting cosmological dynamics happen when the expanding dimensions are still near the 
string scale, and have begun to examine the effects of nontrivial backgrounds 

of fields other than the metric in the low energy supergravity. 

Importantly, several tests have been made of the BV mechanism for determining the 
number^ of macroscopic dimensionality of space. A test of the classical string interaction 
rate |30] agreed with the BV argument that three or fewer large dimensions are allowed. 
However, another early test 3l| of the BV mechanism noted that long strings effectively have 
a width because they can oscillate at a low energy cost, and the effective width invalidates 
the classical BV mechanism. In that event, initial conditions determine the number of 
macroscopic dimensions, which need not be three or fewer. Considering that winding strings 
have a quantum mechanical cross-sectional width also has the same effect, which has recently 
been demonstrated both in 11-dimensional M-theory |32i] and 10- dimensional string theory 
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Ssj l . These tests are all based on the fact that the winding strings will be unable to annihilate 
efficiently if their interaction rate F drops below the Hubble parameter for the expanding 
dimensions. 

In this paper, we will also present a test of the BV mechanism in type II string theory, 
along the lines of 3^.^ However, rather than follow detailed deviations from equilibrium, 
we will present a general argument that winding strings will very rapidly freeze out of a BV 
cosmology. More precisely, we remind the reader that the small-radius, high-energy density 
phase of a string gas is pressureless and demonstrate that any initial conditions which lead 
to an escape of the pressureless phase also lead to freeze out of the strings. Additionally, we 
briefly report numerical studies of a putative high-energy, large-radius phase of a string gas. 
In a second, largely independent section, we calculate explicitly winding string annihilations 
in a BV cosmology, confirming that winding strings fall out of equilibrium immediately. We 
point out that string (or particle) interactions in totally compact spaces generically violate 
conservation of energy, and we use a perturbation theory in the time dependence of the 
background to account for that violation. This technique should be useful in future studies 
of BV cosmology or any cosmology on totally compact spaces. 



II. EQUILIBRIUM COSMOLOGY AND THE DILATON 

In this section we discuss the cosmology of strings in thermal equilibrium in fully compact 
space. In order to do so, we will review the thermodynamics of strings in two regimes; the 
first is high energy with string scale radius, and the second is high energy and large radius. 
In each of these "eras" of string cosmology, we then present solutions to the equations of 
motion. Our discussion of the ffist (string radius) era is largely a review of results presented 



m 



though we will draw some new conclusions based on the behavior of the 
dilaton. 

As is usually the case, we will assume that the spatial dimensions form a rectangular 
torus. For simplicity, we consider two isotropic sets of dimensions, one set which will have 
some initial expansion, and one set which will remain string-scale. For the metric, we take 

ds^ = -de + R'dx'i + 'd^^ . (1) 



^ In fact, [33II appeared as this paper was in the final stages of preparation. 
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where, for the purpose of illsutration, we have assumed that R{t) is the homogeneous scale 
factor for the d expanding dimensions and rest of the spatial dimensions, 9 — d, have scale 
factor R'(t). The coordinate radii are all y/a'. For convenience, we will also define 

i? = , R' = 6" . (2) 

The total volume of the system is given by 

V = {2-K^fR''R! ^-'^ = {iTx^^'fe^'^e'^''-''^ . (3) 

When we discuss specific solutions to the equations of motion, we will generally take z/ = 0. 
Further, we will henceforth take a' = 1. 

Before we start, a few comments are in order regarding our approach to string ther- 
modynamics. First of all, we consider a thermodynamical gas of strings on a classical 
(super)gravity background, as opposed to including the gravitational variables in the ther- 
modynamical treatment The latter approach is mainly important for determining the 
likelihood of different initial conditions; we take a somewhat different tack. Another point 
is that we will use the microcanonical ensemble, which is appropriate for a totally compact 
universe. The microcanonical description of string thermodynamics has been developed in 
HHfly y a. a ..view of *i„g the^^dynan,.., incMng n,o.e .ecent .efe.ence. 
for microcanonical calculations, appears in |4fll |. 



A. Pressureless Phase at String-Size Radii 

First we will discuss the Hagedorn phase of strings when all the spatial dimensions have 
radii near the string scale, i? ~ i?' ~ 1. 

As a preliminary, we should introduce the equations of motion for the supergravity back- 
ground in the presence of a gas of strings. For simplicity, we will ignore the NSNS 3-form 
field strength and all the RR backgrounds, as is customary.^ It is convenient to introduce a 
dimensionally reduced dilaton 

i) = 2^-dfi-{9-d)u, , (4) 



See 



27 



2^ for some progress including the NSNS 3-form and 0, for forms in M-theory. 
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where $ is the dilaton of the lOD string theory. We can now write down the equtaions of 
motion 

- dfi^ - (9 - dy + ^2 = e'^E , (5) 
fl-ijfi = ^e^Pd, (6) 

I'-^O = ^e^Pg-d, (7) 

i>-dfi^-{9- d)u^ = ]^e^E . (8) 

E represents the total energy of the system. The variables P^, Pg-^ are related to the 
pressures pd,P9-d in the respective directions by a volume rescaling, P = pV. Note that 
we will henceforth refer to P as the pressure. (By referring to pressure "in a direction," we 
strictly speaking mean the appropriate diagonal component of the stress-energy tensor.) We 
have assumed that the lOD dilaton $ has no potential in this early stage of cosmology, as the 
dilaton must be free for winding strings to cause the recontraction of expanding dimensions 
. This assumption, of course, does not preclude moduli-stabilizing potentials later in 
cosmology. Additionally, we have assumed that the matter action (which we can think of as 
a thermodynamical free energy) is independent of the dilaton i/); relaxing this assumption 
changes only equation (jH)) [{^. Note that (0) implies that ip is monotonic, and we take it to 
be decreasing, as is usual both to maintain perturbative string theory and prevent runaway 
solutions. Finally, we remind the reader that there are corrections to the equations of motion 
at higher orders in a'; these corrections are negligible in our cosmology as long as all time 
derivatives are smaller than unity. 

Therefore, the important thermodynamical data for our purposes are the total energy 

string scale, 
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and pressure. The appropriate microcanonical thermodynamics for strings in 
totally compact spaces have been discussed in and subsequently in jsj, 35 
The basic thermodynamical relationships in the microcanonical ensemble are 

P_r9S _ ds i_ds 



where S is the entropy. According to js^, the microcanonical density of states at large 



5 



energies (when corrected to account for conservation of charge [35|) is given by^ 

rL+ico dp 



n 



L—ioo 



27ri 



e^-^Z(/3, R, R') , Z 



1 yr / /3o-/3„ y" 



(10) 



On notation: f3o = 2'n{2a'fl'^ = y/s TT^ is the Hagedorn temperature for type II strings, the 
other temperatures are 



2 ^ \Ri 



1/2 



or Po 



1 - 2 ^ 



1/2 



(11) 



for J2{n/Ry,J2{nRy < 2, and gn is the multiplicity of the vectors n with the same 
X](^/-R)^, Z](^-R)^- The ~ indicates equality up to an overall (nearly constant) factor. For 
the radii with 1 < R, R' < \/2, there are four poles, and we get 



l--.r^-^{{m)E) 



2d-l 



-ME 



{2d-l)\ 



{2d-l)\ 



'i7-2dy. 



n-2d 



e-iv'2)E 

[n-2dy. 



{{m)E) 



17- 2d 



2d-l 



(12) 



The poles are at 



A = /3o 



1 1^/2 



2i?2 



/52 = /?0 



R^ 



211/2 



5-1,2 = 2c? 



(13) 



and rji = /3q — Pi. The primed poles are as in (fT^ with R ^ R' and = ~ 2(i. 

Then we obtain the leading order expression for the temperature and the pressures (see 



^ E i=l,2 



j=l,2 



(14) 

(15) 
(16) 



where Ci, C3 are polynomials of t] and C3, C4 are functions of r/'. Note that for large enough 
E, the exponential terms dominate over the polynomial terms, and, as a result, the tem- 
perature remains close to the Hagedorn temperature, T ~ (3q^ and the pressures remain 
vanishingly small, Pd ~ -Pg-d ~ 0. This is precisely the regime of our interest (because 
there should be a macroscopic number of strings), so we will take the pressures to vanish. 



^ The canonical partition function is also calculated in [4lL |42| . 
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Assuming adiabaticity, then, the energy should be a constant (as indeed follows from con- 
servation of the stress tensor). Note that technically we should change the poles once we 
get to R > -\/2; however, at high energies, the pressures are still exponentially suppressed, 
as are the radius dependent corrections to the temperature. Therefore, we can maintain 
the approximation of a pressureless gas that maintains roughly a constant temperature and 
total energy during the first phase of cosmological expansion. 

We should pause to explain why the pressure nearly vanishes. Heuristically, the reason is 
that the winding strings exert a negative pressure (like cosmic strings), and, with all radii 
near the string scale, the population of winding modes and momentum modes are nearly 
equal (because their masses are very similar). In fact, T-duality implies vanishing pressure. 
At the self-dual radius, the pressures P, P in the original and T-dual variables should be 
equal. However, because /i, — > —fi, —v under the associated T-duality, P = —P, so the 
pressure vanishes. Therefore, even more general gases of strings and D-branes must have a 
pressureless phase near 

In this pressureless phase, the equations of motion are relatively easy to solve, and in fact 
are a simple generalization of results given in 11^. Isi | : 



^t' + Bt 4 
4 

/io H log 

a 



I ^^-(i 1 
uo H log 

a 



- dAl - (9 - d)Al_, 
Eo 

'{Eot + 2B -2a){B + a)' 
{Eot + 2B + 2a){B - a) 
{Eot + 2B -2a){B + a) 
{Eot + 2B + 2a){B-a) 



a 



dAl + (9 - c/)A 



9-d ) 



(17) 
(18) 

(19) 
(20) 



where 



/ioe 



-i>o 



A 



9-d 



-V'o 



B 



-i'oe' 



-V'O 



(21) 



A useful intermediate result, which follows directly from equations ()fif7j) is that both /i, z> cx 
e"^. In fact, this result is true in the pressureless phase even if (jHJ is modified. For simplicity, 
we will consider only solutions in which z/(t) = 0, so Ag^^i = 0. We henceforth drop the 
subscript on A^. See figure [T] for plots of the scale factor and dilaton in the distinct cases 
of (i = 3, 9 expanding dimensions; note that the time evolution of ijj is not that sensitive to 
dimension. 
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(a) fi vs t 



(b) ip vs t 



FIG. 1: Cosmological evolution of the scale factor and dimensionally reduced dilaton for fiQ = 0.7, 
■00 = —4, and £"0 = 100 for 3 (solid curve, diamonds) and 9 (dashed curve, stars) expanding 
dimensions. 



As discussed in [l^, the radius R = asymptotes to the value 



Rr, 



(22) 



B - VdA 

Because the radius asymptotes to a finite value, the only hope that we have of "decompact- 
ification" of the fi directions is for the string gas to leave the pressureless phase. As we have 
noted in our discussion of the thermodynamics, the gas of strings will leave the pressureless 
phase when the universe is no longer near the string radius, when > R for some radius R 
which divides "large" from "small." We take ^ ^ 3 for numerical purposes. In particular, 
for the asymptotic radius i?oo > -R, we find that 

2 

- 1 



{Re 



MO 



+ 1 



{Re- 



MO 



1 



(23) 



As far as we know, this is the first discussion of such a limit on the initial dilaton. Note 
that at large energies, the dilaton is considerably suppressed. 

Finally, we should check that our initial conditions avoid a Jeans instability resulting 
in a black hole, as was discussed for string thermodynamics in To avoid the Jeans 

instability, we require 



R' < 



1 



(24) 
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For all radii initially at (j24j) is 

e* < ^ . (25) 

The constraint is somewhat less restrictive than if we keep /io < 1 as required for 
the use of the low energy equations of motion. 



B. Hagedorn Strings at Large Radii 

Now we would like to discuss the cosmological evolution associated with a phase of strings 
at high energies and large radius, which we take to start at R > R. Although we will find 
later that the initial, pressureless phase cannot remain in equilibrium in order to lead to 
this phase, presumably strings could equilibrate after the universe reaches a large radius. 
Therefore, these cosmological solutions may be useful in future studies. We know of no other 
reference that has considered the cosmological evolution of this phase of strings, as other 
works have concentrated on lower energy gases of strings, which are dominated by radiation 



see 



i4 



For this later era 



give 



(26) 



for the multiple string density of states. We include the prefactor of because momentum 
and winding in each of the nine spatial directions are conserved, as discussed in 
Here A is a constant (which we take to vanish), and a is constant with radius but depends 
on the number of large dimensions: 



_2d/2 



27r'^/^ 
r{d/2) Jo 



dxx'^ ^ In 



1 + ^2 



(1-.^) 



1/2 



(27) 



The entropy is therefore 



S = \xi(3q-Q\yiE + (3oE + {\ + a)R'^ 



(28) 



The temperature and pressure are therefore 



1 

T 
P 



PoE-9 
E 



(A + a)de^"^ . 



(29) 
(30) 
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(a) fi vs t 



(b) ip \s t 



FIG. 2: Evolution of the scale factor and dilaton in the large-radius Hagedorn phase for d = 3 
expanding dimensions. The initial conditions are given by the pressureless phase evolution of the 
initial conditions given in figure ^ 

Working with the usual assumption of adiabaticity, the equations of motion are now too 
difficult to solve analytically. The main stumbling block is that, holding the entropy ()28|) 
fixed, the relation between the energy and radius is transcendent. However, it is straight- 
forward to solve the equations of motion numerically; we found it useful to rewrite 



and solve numerically for E function of time. 

We include plots of the scale factor and dilaton for this large-radius phase in c? = 3 
expanding dimensions in figured For initial conditions, we take the values of E, fi, ip, and 



might be more interesting, as indicated below). In the pressureless phase, we took the same 
initial conditions as listed in figure ^ Now the expanding dimensions can resume a rapid 
expansion, while the dilaton continues its monotonic decrease. When the energy density has 
decreased enough due to cosmological redshifting, this phase should match onto a radiation 
dominated phase (assuming thermal equilibrium). 




(31) 



ip given by evolving the pressureless phase to e'^ = -R = 3 (although other initial conditions 
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C. Estimating Interaction Rates from Cosmology 



Now we will consider the importance of (j^^ in determining the interaction rate of strings. 
As a particular interaction of interest for the BV mechanism, we can focus on the annihilation 
of winding and anti-winding string pairs, although the key issue is whether the strings can 
remain in equilibrium. We first give a rough estimate of the interaction rate. 

Up to gamma-function factors, the amplitude for any 2 string to 2 string process (with all 
dimensions near the string radius) is e"^, as we will discuss in more detail below. Therefore, 
we write the interaction rate roughly as 



where A/" is the number of strings with which a given single string can interact, and the 
sums are over the outgoing spin and momentum states. At this level of discussion, the sum 
over spins just gives a constant factor of order 100. Also, for this simple estimate, we can 
approximate the sum over outgoing momentum states as a constant factor; for dimensions 
near the string radius, corresponding to the pressureless phase, this momentum sum factor 
is order unity. 

Since the classical BV mechanism distinguishes the annihilation rates of winding strings 
in (i < 3 and d > 3 large dimensions, we comment here on the effect of varying d on ()32|1 . the 
quantum mechanical annihilation rate. The dimensionality enters into F through the time 
evolution of ip (and also somewhat in the constraint on its initial value (I2S1)); the higher d, 
the faster ijj decreases. This effect tends to enact the BV mechanism, as strings interact less 
rapidly with more growing dimensions. The other factor through which d affects F is the 
sum over outgoing momenta. Quite simply, for a given energy, the phase space of outgoing 
momenta is larger in a greater number of large dimensions. This tends to counterbalance 
the variation of the dilaton, but it is a small effect when all dimensions are near the string 
radius. 

The number of strings available for scattering is determined by thermodynamics. For the 
pressureless phase, the number of strings of charge q and energy between e and e + cie in a 
gas of strings of total energy E is 



F^A^5:5:e^^ 



(32) 



X{e,q,E) 



Po) e \e{E-e) 



) 



exp 



Ae{E - e) 



E 




(33) 
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Here the matrix A is 



A = -A^diag(l//?, ■ ■ -, 1/R', ■■■,/?,■■■,/?',■■■). (34) 



(4vr) 



d 9-d d 9-d 

This takes into account the contribution of winding and momentum, as well as oscillator 
modes, to the energy of the string. The total number of strings in the gas, given by summing 
over charge and integrating over e is A/" = / deJ2q-^{^, E) ^ Inii^. Actually, restricting to 
strings of a certain charge or range of energies gives a considerably smaller number, so we 
can take A/" ^ In as an upper bound in (j32j) . 

If we substitute in the bound (j^ . our estimate of the interaction rate becomes 



r< 100(^1^^ InEo (35) 

during the entire pressureless phase. We have used the facts that the total energy in the 
string gas is constant during the pressureless phase and that the dilaton is monotonic de- 
creasing. Even for Eq ^ 100, which only leads to about 4 strings in equilibrium, we have 
r < /io- For larger energies, the interaction rate will only decrease. 

Therefore, it seems inevitable that the string gas will fall out of equilibrium almost 
immediately as the universe begins expanding in the pressureless phase. This means that, 
for any d, winding strings will freeze out and come to dominate the energy density of the 
compact universe and cause its recollapse. Contrast this case to the original BV mechanism, 
which claims that winding strings can annihilate effeciently for d < 3.^ The reader might 
wonder if the string gas could come into equilibrium as /i decreases in time; however, if we 
assume pressureless evolution, ft oc while F oc e^'^ decreases more rapidly. 

As a result, any study of the standard BV scenario should focus on the non-equilibrium 
thermodynamica and cosmology of string gases, as partially discussed in 0] (these refer- 
ences included the non-equilibrium Boltzmann equation for strings but used the equilibrium 
cosmological evolution). There are a few possibilities to avoid this conclusion; however, they 
all modify the usual BV scenario. A simple idea would be to consider lower energy string 
gases (ie, /? ^ Pq), as in j20i], but winding states would be populated only infrequently, 
so it might be difficult for winding strings to stabilize the radii of any of the dimensions. 



^ Of course, in more complicated compactifications, such as orbifolds, winding strings can decay _oiit_ of 
equilibrium. Then the question is whether they decay too rapidly to produce the BV mechanism 
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One obvious possibility is to take fiQ ^ 1, which would require understanding the full time- 
dependent worldsheet CFT. Unfortunately, this is beyond our capabilities, and in fact such 
a CFT may not have a straightforward geometric interpretation. Another possibility would 
be to consider the monotonically increasing solution for the dilaton, which could lead to a 
strong gravity regime and large time derivatives. Alternately, heterotic strings have qualita- 
tively different thermodynamics due to additional low-energy states near the string radius, 
so they may evade (j23|) . However, at high enough energies, even heterotic compactifications 
should have a phase similar to the pressureless phase described above. 

Another attractive option is the introduction of winding branes in the gas of strings, 
as has been widely considered 5, 0, Q, 0]. Gases of D-branes might avoid our constraints 
because the branes couple to the dilaton and therefore modify equation (jH)). However, the 
other equations of^motion remain the same js], and, because brane gases are indeed T- 
duality invariant ^fij, the early phase of cosmology should still be pressureless. Therefore, 
/i oc e'^ as before. As long as ip decreases quickly enough, ^ will still asymptote to a 
finite value, leading to a constraint on the initial value of the dilaton, as in (j^Hj) . In a 
perturbative model, in fact, we expect D-brane states to be populated very little, so branes 
will not have much effect on the time evolution of ip. Therefore, it seems reasonable that 



will noia approximatel y ev en 

jnl Q 3, HH Q 21, 22, 



hold a ppr oximatel y ev en i 



branes are included in the thermal gas of strings. See 
261 1 for some work relevant to time evolution in brane 
gases. Similarly, allowing the shape moduli of the spatial torus and the various form fields 
of supergravity to be nontrivial may alter the qualitative behavior of the dilaton, such as 
)y providing a potential (for a review of flux compactifications with stabilized dilaton, see 
4J| ) . Some work including other deg rees of freedom in the context of string and brane gas 



cosmology appears in 



271 . l28| |. We must caution that including a potential for the 



dilaton appears to prevent radius stabilization by winding strings, however 24 1. 

As a final alternative, we propose a stochastic version of the BV scenario. In this case, we 
could imagine that the dilaton is large enough to maintain thermal equilibrium of strings (so 
that the winding modes can annihilate), but the expanding dimensions asymptote to a small 
radius. Then, if some of these dimensions were to start expanding again due to some thermal 
or quantum fluctuation, /xq > 0, so the final radius can grow. Eventually, though a 
random walk, some dimensions could grow large while maintaining thermal equilibrium in 
the string gas. There are some difficulties with this proposal, however. One is understanding 
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the quantum mechanics and thermodynamics of the supergravity background. Another is 
that the dilaton will still decrease monotonically in time (except perhaps for thermal and 
quantum fluctuations), so the string gas might fall out of equilibrium in any event. 

All of the above proposals are candidates to allow the string gas to remain in thermal 
equilibrium and the winding modes to annihilate. However, it is far from obvious that any 
of them discriminate among the number of expanding dimensions d. That is, in any of the 
above proposals, the values of d which allow strings to remain in thermal equilibrium until 
all the winding strings annihilate (statistically speaking) seem likely to depend on initial 
conditions rather than d itself. Perhaps the most encouraging of our proposals in that 
regard is the inclusion of other form fields, which have varying rank, which could influence 
dimensionality to some extent. 



III. INTERACTIONS IN TOTALLY COMPACT SPACES 

In this section, we will calculate the annihilation rate for a winding and anti-winding 
string pair in order to confirm that the strings cannot stay in equilibrium in the pressureless 
phase. In so doing, we point out a problem with the naive use of string perturbation theory 
in the BV scenario and point out a resolution that could be useful in future tests of the BV 
mechanism. 

As we have discussed, a key feature of the string-gas and brane-gas cosmologies is that all 
the spatial dimensions are compact; in fact, much of the interesting dynamics that determines 
how many dimensions can expand to macroscopic size occurs when all the dimensions are 
within an order of magnitude of the string scale (loitering, for example, occurs at roughly 
2 to 3 times the string length 2^ 13] )■ Because all the dimensions are compact and in 
fact at a small scale, momentum must be considered as being quantized. When the large 
dimensions become much larger than string scale, we are justified in ignoring the momentum 
quantization because the energy spacing will be much smaller than other energy scales. 
However, since we are interested in string scale dimensions, the energy spacing of momentum 
modes is nearly as large as the string mass or even winding modes, we have to quantize 
momenta, even those in the "large" dimensions. 

Quantization of all the momenta puts interesting constraints on the interactions of strings, 
and, in particular, at generic compactification radii, interactions cannot conserve energy. In 
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other words, superstrings on a time-independent (or presumably any other compact 
space) cannot interact. To see this in more detail, consider some simple string amplitudes: 

1. Supergravity mode interaction in which two strings with momentum in the direction 
scatter into the direction. In the simplest case, the strings have no other momenta, 
so the initial strings have rii, —rii momentum numbers and the final strings have 
^2, —712- Conservation of energy would then require 2ni/Ri — 2n2/ Ri, which is clearly 
impossible at generic radius since the e Z. 

2. Oscillator mode decay to momentum modes. If the final state strings have momenta 
in one direction only, energy conservation would require N — 2n/ R, which is violated 
at nonintegral radius. 

3. Winding mode annihilation, which is our main focus. In the simplest case, there are 
two strings with opposite winding number w, —w in the x^ direction, which inter- 
act to produce, for example, momentum numbers n, —n in the x"^ direction. Energy 
conservation would be 2wRi — 2n/i?2, which is again clearly impossible for generic 
radii. 

In fact, even at radii (for example, R — 1) where the simplest case interactions listed above 
are allowed energetically, there are string states with multiple nonzero quantum numbers 

that cannot interact without violating energy conservation. Essentially, the problem is that 
at any given radius only a small subset of string interactions might be allowed by conservation 
of energy. 

There are several possible resolutions to this difficulty. One is that the radii are quantum 
mechanical variables, so that there is some probability of being at a radius where any given 
interaction is allowed kinematically. Second, the more energetic states might have a width, 
or uncertainty in their energy, due precisely to the interactions in question. The logic is 
somewhat circular, though, as the width allows the interactions which give rise to it. A 
final solution to the dilemma is to remember that energy is not strictly conserved in an 
expanding universe. For a scale factor given hj R — e^, we expect that there will be energy 
nonconservation of order fi in each interaction. Since the smallest energy quantum is 1/R, 
we expect that all interactions would be allowed when the radius of the universe is longer 
than the Hubble time. We will focus on this third solution in this section, as it is in some 
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ways the most straightforward to understand, and we discuss a formahsm for deaUng with 
compactification radii quantum mechanically in the appendix. In all of these approaches, 
we would assume that, although individual interactions would not conserve energy, the total 
thermodynamical energy would be conserved (that is, the nonconservation of energy would 
average to zero over many interactions). 

A. Time Dependence and Energy Nonconservation 

Because we have been careful to ensure that all time derivatives in the supergravity 
background are smaller than string scale, we can use a very simple approach to include 
the expansion of the universe in the string interactions. Namely, rather than incorporate 
time dependence directly into string perturbation theory, we develop an effective quantum 
mechanics based on string perturbation theory and then consider the time dependence of the 
effective quantum mechanics. We will focus on winding mode annihilation to momentum 
modes, WW NN, but our approach can easily be generalized to any interaction. 

In its simplest form, we can take a two-state quantum mechanics, with the winding/anti- 
winding string pair (PFiy) as one state and the momentum string pair (NN) as the other. 
We derive the Hamiltonian as follows. The diagonal, mass part comes directly from the 
perturbative string spectrum on the fully compact space, and we do not need to discuss it 
further. Including an interaction piece, the Hamiltonian is 



assuming that the momentum and winding vectors n and tv are completely in the expanding 
directions (the factors of 2 in the diagonal elements are because each state represents a 
pair of strings). To find Vi, we compare to the corresponding amplitude 27c6{Ef — Ei)A for 
WW — s> NN in string perturbation theory (which we discuss in more detail in section UlI B|) . 
Standard string perturbation theory assumes a time independent background and allows an 
infinite time for the strings to interact. In quantum mechanical perturbation theory in this 
case, we find the same amplitude 



2\w\e>' 



Vi 



H 



(36) 



Vi 2\n\e~^' 





(37) 
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for a time independent Vj. Taking into account the relativistic normalization of states in 
string perturbation theory, 

Vi = - I / , (38) 

where the Ei are the energies of the incoming and outgoing strings. Extracting the depen- 
dence on supergravity fields, Vj = Ce^, where 

C = '"^-^ (39) 

^(2Ei) (2^2) (2^3) (2^4) 



is a constant factor. One advantage of this effective quantum mechanics approach is that the 
string amplitude automatically includes single string intermediate states, as we will discuss 
in more detail below. To include the evolution of the universe at any given time to, we can 
expand H to first order in time derivatives, H = Hq + Ht. Then we can treat the first order 
term as a perturbation around the zeroth order term. 

We should therefore first find the eigenstates of the time independent Hamiltonian, in- 
cluding the interactions Vj. In this way, we avoid confusing state oscillation of the purely 
time independent quantum mechanics in the annihilation rate, which should be due purely 
to time dependence in the Hamiltonian. Using standard formulae from quantum mechanics 
(see, for example, eigenstates are given by 

\w,-w) = \w, + \n, -nf 

w n 

|n, -n) = \n, + — -wf , (40) 

n w 

where the energy eigenvalues are just the unperturbed ones E^, E^. We have treated Vi as 
a perturbation because we are working in a regime in which string perturbation theory is 
valid (ip <0). 

n 

The lowest order transition amplitude for winding annihilation is (again, see |45[, for 
example) 

c(^) = -i( Ce^i; + -2Ce^wfxe^' 

\ UJ 



UJ 



f-w- (41) 

lUJ [ZUJ)^ 



with 

2n 

u = Ef - Ei = — - 2wR . (42) 
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Here, n = \n\ and w = \w\. 

Next we note that fit, tpt < 1 in order to ensure that we can treat the time dependence 
as a perturbation. In fact, this is consistent with assuming t ~ 1 for the interaction time 
of strings, as the supergravity approximation requires that all time derivatives be less than 
unity. Then we can take ut to be small, so the probability is therefore 



\Ce^\' 



ip"^ wfie^ip fine ^ip 
4 u uj 



(43) 



The total annihilation rate for a winding string {WW NN), as discussed in section 
III (^1 following equation (jH^ . is 



w,n) 



(44) 



with notation as in III CI Since we will ignore string polarizations for simplicity, the sum over 
outgoing spins reduces to an overall constant. In fact, conservation of angular momentum 
correlates the two outgoing spins, so the overall factor will be the number of massless super- 
string spins, 256, times a symmetry factor of 1/2. The sum over outgoing momenta should 
be within the (larger) expanding dimensions, as that is the type of final state that inter- 
ests us. Before we present the results of a detailed computation, we will discuss the string 
perturbation theory amplitude A that we used to define the effective quantum mechanics. 
One note: the reader may wonder why we did not take advantage of the optical theorem 
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47|. The simple 



to calculate the total interaction rate (up to the factor of A/"), as in 
reason is that we have no simple way to incorporate the time dependence of the background 
into the appropriate calculation in string perturbation theory, unlike the direct calculation 
that we present. 



B. Bosonic String Winding Annihilation Amplitude 

As we discussed above, we need the string perturbation theory interaction amplitude. 
We carry out a somewhat rough analysis of the interaction, ignoring the contribution to 
the amplitude from string polarizations. This approximation allows us to replace the sum 
over outgoing spins by an overall factor (and the average over incoming spins by unity) 
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in the interaction rate (j44j) . Therefore, we use the bosonic string amphtude for nonpolar- 
ized strings (the tachyon mode if the winding and compact momenta vanish) but with the 
appropriate mass-shell for whatever polarized superstring modes we are considering. An 



analogous superstring calculation for tachyonic winding strings was discussed in j47| using 
the optical theorem; we calculate the amplitude directly at tree level for use in the effective 
quantum mechanics as discussed above. Our presentation is a generalization of the answer 
to a textbook problem from 3|. 

In calculating the annihilation amplitude, we find it convenient to work with a unit metric 
and coordinate radii equal to the proper radii of the spatial dimensions, but the amplitude 
is diffeomorphism invariant and can therefore be used directly with the conventions in the 
text. The vertex operators of the four interacting strings are all 

: exp [iki ■ Xl + ik^ ■ Xr] : , (45) 



where the closed string coupling 

^c=^(27r)^/W, (46) 

is the lOD dilaton, and V is the total volume of the spatial dimensions (the prod- 
uct of 2ttR over all dimensions). The amplitude is just the expectation value of four 
vertex operators, and it has a prefactor (besides the vertex operator normalization) of 
ii^STr/ g'^a'){2n)6{J2 k^)VSnSw Here S(J2k^) conserves energy and 6n,w conserve compact 
momentum and winding (ie, these are products of Kronecker delta symbols). There is also 
an overall sign (the cocycle) which we will ignore as we have only a single amplitude. Thus, 
the amplitude becomes 

A = z2{27i)-h^{2n)6'^~\k)6J^F{kL, kn) , (47) 

with ip as before. Here F is some function of the momenta that we need to determine, and 
we are again setting a' = 1. 

As it turns out, we are interested in annihilation of winding, so the two outgoing strings 
will have zero winding, and the two incoming strings will have opposite winding w = wi = 
—W2- In addition, we will take the simplest case in which the winding strings have no 
momentum. The compact momenta are 



^L3 — — —kiA — —km 



khi = -km = -kL2 = kR2 = (wR) (48) 
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where {n/R) is the vector with components rii/Ri {Ri is the physical radius in the ith 
direction) and similarly for (wR]. Additionally, kj represents the incoming energy of the 
Ith string, which is negative for the outgoing momentum strings, while Ej is the physical 
energy of the Ith string, which is always positive. 
In this case, the momentum dependence becomes 



/ 



d'zlzf^^^ll - zf-^^ ( - ) ( ^ ) (49) 



.zj \l — z, 

integrated over the sphere. Note that £'3^4 appear asymmetrically due to the way in which we 
have located the vertex operators on the sphere; however, the end result will be symmetric, 
as it must due to SL(2, C invariance on the worldsheet. To do this integral, we can integrate 

by parts n-w on each of z. Then we some prefactors along with the usual integration needed 
for the Virasoro- Shapiro amplitude. As it turns out, the F functions from the integral can 
absorb the prefactors, giving us 

F — F1F2F3 



Fi 



r(-i(EiE4 + n-w)) r(-i(nw + 



n ■ w 



r(l + l(E2E4 + n-w)) r(l + l{nw + n-w)) 
r [-1{E2E^ - n ■ w)) ~ r[-l{nw-n-w)) 
p ^ r{-l-\{E, + E2)E,) ^ r (-1 - nw) 
' t{2 + \{E, + E2)E,) T{2 + nw) ■ ^ ' 

The second expression for the Fj arise from using the mass-shell condition for the strings, 
which should be of the form (in the general case) 



-i>2 



£;2 = 4Ar+g) +{^\ (51) 

is an integer that describes the oscillator excitation of the string. For bosonic string 
tachyon vertex operators, such as we are using, we should take A" = — 1; however, we want 
to simulate superstrings at the massless level, so we take A" = 0. We have therefore taken 
Ei^2 = wR and £^3^4 = n/R. 

There is one additional subtlety to discuss. By unitarity, the string amplitude contains 
poles corresponding to single string intermediate states. (This fact is actually very useful for 
our quantum mechanics because it means we do not need to use second order perturbation 
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theory to include all interactions at the same order in e'^. All the tree-level annihilation 
processes are included in this single sphere amplitude.) These poles appear through the 
gamma function in the numerator of F3, which has a pole whenever nw is an integer. Indeed, 
if we rewrite nw in terms of the center-of-mass energy (Mandelstam variable) nw = s/4 = 
{El + E2Y = {E^ + E^y, then the poles are at s = —4, 0, 4, 8, ■ ■ ■, just the bosonic closed 
string oscillator spectrum. 

However, since we are interested in computing the actual interaction rates, we must 
include the widths of the poles. Although we do not know the appropriate formulation of 
resonances in the full time-dependent formalism, we can, at our level of approximation, use 
the traditional Breit-Wigner resonance, replacing s — > s + imTy^ at the pole given by s = m^, 
where is the width of the resonance (see a quantum field theory text, such as To 
encompass all the poles at once, we use s s + iy/sTw, which is correct at each pole. To 
calculate the width, we ignore string polarizations again, so the width becomes 

r» = i28i:5^;(4'*«(v^-2"/fl) . (52) 

The delta function conserves energy, and ^/s is the energy of the intermediate oscillator 
resonance. Due to momentum conservation, the sum over outgoing momenta is a single sum 
in the d expanding directions, which we approximate by an integral: 



Therefore, the width is 



leaving us 



GAR'^e^ ( - I ^ ^ (54) 

U) r((rf-l)/2) ' ^^^^ 



T{-l-nw -i^/muT^/2) 

^3 = wjT^. ^ , (55) 

r (2 + nw) 



with s replaced by Anw in F^. 



Results of Calculation 



We now briefly present the results of calculating the total annihilation rate ()44p in the 
pressureless phase of the cosmology. For the purposes of the calculation, we first take 
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FIG. 3: The interaction rate for strings in the pressureless phase for d = 3 expanding dimensions. 
The initial conditions are those given in figure ^ The lower, solid curve has M calculated for 
w = —1, and the dashed curve takes M = lnE. 

Af = 7V(e = wR, q = w, E)Ae, where we take Ae = 1 in string units. Thus, we are 
calculating the annihilation rate purely for strings with w units of winding and no extra 
oscillators. Then, as an upper bound on the annihilation rate, we take Af = InE, which 
is the total number of strings in the thermodynamical gas. We see that, in either case, 
the interaction rate is small compared to the Hubble constant /i, so the gas must fall out 
of equiUbrium. We do the momenum sum Y.n by brute force, requiring at all times that 
the frequency uj = 2n/R — 2wR < 1 (since larger values are suppressed). For our example 
calculation, we take w = (1, 0, ■ ■ ■ 0), and i? < 3 in the first era, so we take each component 
of n between —10 and 10 for the sum over momenta. The logarithm of the interaction rate 
is plotted in figure El The initial conditions are the same as those used in plotting figure d 
Note that the gamma functions give an additional suppression of the interaction rate. 

The interaction rate is very small compared to the Hubble parameter (which is order unity 
as an initial condition). We have carried out this calculation also for d = 2,4 expanding 
dimensions and found no qualitative difference. There is one interesting feature of these 
interaction rates: a slight increase at later times. This temporary increase is due to an 
increase in phase space available to the outgoing modes as some of the dimensions expand. 
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IV. SUMMARY 



In summary, we have made two arguments regarding the usual Brandenberger-Vafa pro- 
posal for the cosmology of string gases on totally compact spaces. 

In the first part of the paper, we analyzed two possible states of a string gas at high 
energies, reviewing their thermodynamics and solving for the cosmological evolution of the 
universe in those phases. The later phase had not been previously studied in the literature. 
From the time evolution of the dilaton and scale factor, we were able to give a heuristic 
argument that there are not initial conditions in the small-radius, pressureless phase that 
both grow into the large-radius phase and maintain thermal equilibrium of the string gas. 
We also argued that more general brane gases would likely have a similar constraint and 
gave several possible scenarios that could avoid this conclusion, all of which modify the 
usual Brandenberger-Vafa scenario. One key point is that our conclusion that strings fall 
out of thermal equilibrium is largely independent of the number of expanding dimensions. 
Therefore, it would seem like any possible modification to the Brandenberger-Vafa scenario 
would not select a favored number of large dimensions, as was originally proposed in Q| . 

In the second part of this paper, we calculated the approximate annihilation rate for 
winding strings in different numbers of expanding dimensions. This calculation confirmed 
that the strings cannot remain in thermal equilibrium in the small-radius phase of the string 
gas. In that section of the paper, we also argued that quantization of string momenta is 
important during the small-radius phase. Specifically, due to quantization of string wind- 
ing and momenta, string interactions will generically violate conservation of energy (and 
would thus seemingly be forbidden). We demonstrated that the time dependence of the su- 
pergravity background allows string interactions even at small radii, using time dependent 
perturbation theory in an effective quantum mechanics. We believe this approach will be 
useful in future studies of small-radius Brandenberger-Vafa cosmology. 

In conclusion, then, we have argued that it seems unlikely that string gas cosmology 
in totally compact spaces provides an explanation for the number of macroscopic dimen- 
sions, especially when our work is taken in concert with 0, 2^ |^. Nonetheless, the 
Brandenberger-Vafa proposal is a very natural initial condition for cosmology from the 
point of view of string theory, so it is worthy of continued study. We have in fact given some 
directions which we believe would be interesting for future work. 
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APPENDIX A: QUANTUM RADII 

In this appendix, we would like to sketch another approach to resolving issues of energy 
conservation in completely compact spaces based on elementary considerations of quantum 
gravity. In quantum gravity, we should really consider the metric to be some quantum 
variable, so that the radii of the compact dimensions should not be fixed to the classical 
cosmological trajectories. Rather, there should be a wave function ^(//, i/) giving the prob- 
ability density for different values oi /ijU, which should be peaked at the classical trajectory. 
In fact, in the perturbative string regime, we expect that these are sharp peaks because the 
Planck mass is very large (even compared to the string scale). (We expect that ^ is most 
simply written in terms of /i, u because those variables have canonical kinetic terms.) 

Then the total state of the system (at fixed time) should be described as | {matter}, /i, i/), 
where the matter part describes whatever string might be currently in the universe. For 
simplicity, all but two of the strings will be taken to be spectators (as is usual in perturbation 
theory), and they will give only a factor of unity. Therefore, the total amplitude A — 
({matter}, /i, i/|{matter}, /i, v) should be 

J diidu\^{lj,u)\^A6(j2Eil^^'^)) ■ (Al) 

Here, ^ is a reduced string amplitude; just the amplitude from string perturbation theory 
with the energy conserving delta function extracted. We will take the wavefunction to 
factorize on /i, u. To simplify notation, we will henceforth let E be the sum over energies 
(with outgoing energies taken to be negative). As an aside, we are ignoring any complications 
to the measure caused by our constraining sets of dimensions to have the same radii. 

Then suppose we have the simple case described above for winding mode annihilation. If 
both the winding and momentum directions have the same radius (ie, both are among the 
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large dimensions represented by /i), we can evaluate the integral easily to get 



1 



- In — 

V2 w 



) 



2 



A 



A . 



(A2) 



Here we have integrated over the v dependence of the wavefunction. Other examples proceed 
in the same way, with a slight complication if the energy depends on both and v. In that 
case, we will find /i as a function of v by integrating over /i first. Then generally we will 
have 



where /i is the value of that satisfies conservation of energy as a function of u. If the 
wavefunction is sharply peaked in then 



where \1/ now depends only on and everything else (including Jj) are evaluated at the peak 
(classical) value of v. 

The interaction probability should be given by the normalized square of the the amplitude, 
as usual. To get the rate, however, we must divide the probability by amount of time over 
which the interaction has been allowed to take place. In perturbative string theory, as in 
perturbative field theory, we have essentially assumed that strings have been allowed to 
interact over some long period of time. Usually, the square of the amplitude has two energy 
conserving delta functions, and one of those delta functions, evaluated at zero, becomes the 
interaction time. In our analysis, the radial wavefunction replaces the energy conserving 
delta function and should also represent the interaction time. 

We can understand this fact by considering the classical limit of our amplitude. If we 
turn off quantum gravity (perhaps by sending the dimensionally reduced dilaton if) — oo), 
then we should have |\l/(/i,z/)p = 6{fi — fi)6{i' — 0), where /t, z> are the classical values. Then 
we obviously get the usual amplitude. However, we can also integrate over /i, u, so the 
amplitude is as in ()A4jl . Then we note that 



Therefore, when we square the amplitude, we have two factors of \'$\'^/{dE/dfi), one of 
which should be evaluated at the classical values fi, i) and divided out to give the rate. 




(A3) 




(A4) 




(A5) 
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Once we average over incoming polarizations and sum over outgoing states, we end up 
with a rate 

Here the notation is as in equation ()44|) . and we have speciahzed to the 2-string annihilation 
or scattering rate. 

What is missing from this prescription is a wavefunction for the radii, as well as its 
dynamics. One particularly important question is how each string interaction affects the 
wavefunction; that is, does the string interaction "collapse" the wavefunction or is there 
some more complicated entanglement or decoherence process? Another point concerns the 
quantum mechanical evolution of the wavefunction. For example, if we choose to model the 
wavefunction as a gaussian, how does the width of the gaussian evolve with time? 
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